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پایش و پیش‌بینی خشکسالی در استان سیستان و بلوچستان 


چکیده 


خشکسالی کشاورزی در ایران به عنوان یک فاجعه‌ی طبیعی برای تمام مردم» به ویژه 
کشاورزان» شناخته شده است. هرچند از وقوع خشکسالی نمی‌توان جل و گیری کرد اما اگر 
ماهیّت و خصوصیّات آن مورد مطالعه و بررسی قرا رگیرد» می‌توان نسبت به پیش‌بینی آن امیدوار 
بود و با آمادگی از اثرات زیان‌آور آن کاست. در این نوشتار استان سیستان و بلوچستان با توجه به 
داشتن شرایط خاص اقلیمی از نظر پدیده‌ی خشکسالی» مورد بررسی قرارگرفته است. دراین 
پژوهش» پایش و پیش‌بینی خشکسالی استان سیستان و بلوچستان, با استفاده ازالگویی که بر اساس 
شاخص 1و درصد از نرمال و زنجیره‌ی ما رکف تهیه شده .انجام گرفته است. نتایج حاصل از 
تحقیق نشان می‌دهد که منطقه‌ی مورد مطالعه نسبت به خشکسالی‌های ميان مدت و کوتاه مات 
حسّاس است و بیشتر خشکسالی‌هایی که در این منطقه اتفاق افتاده است» خشکسالی ملایم و 
متوسط است و احتمال وقوع پدیده‌ی خشکسالی در مناطق م رکزی استان پیشتر است. 
کلیدواژه‌ها: خشکسالی؛ ۳1 شاخص درصد از نرمال زنجیره مار کف.سیستان و بلوچستان 


۱. درآمد 


یکی از مظاهر طبیعی متأثر از وضعیّت آب و هواه خشکسالی است. به گونه‌ای که کاهش یا افزایش 
تأثیریک با چند عامل از عوامل اقلیمی موجب پیدایش این پدیده می‌شود. خشکسالی کشاورزی در ایران 


به عنوان یک فاجعه‌ی طبیعی برای تمام مردم» به ویژه کشاورزان» شناخته شده است (ارشادی»۸۶. 


۲ مجلهی جغرافاوتوسع‌ی نایا شمارمی‌بازدهم ‏ 

کارهای اولیه در خصوص پایش عمدت با تحلیل فراوانی منطقه‌ای خشکسالی در مقیاس کوچک آغاز 
شد(۷۷[0016,1966) و سپس توسط سایر محققین برای ایالاتی از آمریکاء که بیشتر در معرض 
خشکسالی بودند» مانند ایالات غرب و جنوب غربی دنبال شد(0120100,1986:500,1980)). هیز و 
همکاران (21.,1998 64 1165 ) خشکسالی گسترده سال ۱۹۹۶ در سطح آمریکا را با استفاده از شاخص 
1 پایش نمودند. نقشه‌ها گونه‌ای تهیه شد که وضعیّت عمومی هر کدام از ایالات به طور یکجا تعیین و 
نهایتً قشه‌ی کل کشور تهیه شد. رضیتی و همکاران(۱۳۸۲) نیز پدیده‌ی خشکسالی را در استان‌های 
اصفهان و یزد با استفاده از شاخص 5۲1 و برای دوره‌های زمانی ۳ ۶ ۱۲ و ۲۴ ماهه برای سال‌های ۱۹۹۸ و 
۹ مورد مطالعه قرار دادند. 

پیش‌بینی خشکسالی پر از مشکلات و مسائل خاص است. به طور کلی» روش‌هایی را که توسط آنها 
می‌توان پیش‌بینی‌هایی برای رخدادهای آتی خشکسالی نمود» روش‌های آماری» آماری-فیزیکی» 
سنجش از دور و سینوپتیکی هستند. آنچه در این تحقیق مورد استفاده قرار گرفته‌است. کاربرد روش 
آماری زنجیره مار کف است 

گابریل و نیومن(۱۹۶۲)" اولین کسانی بودند که از یک زنجیره ما رکف مرتبه‌ی اول برای تعیین احتمال 
وقوع بارندگی روزانه در تلاویو استفاده کردند. اوچلا و همکاران(۲۰۰۳) ۰ از اگوی زنجیره مار کف 
مرتبه‌ی اوّل برای شبیه سازی و پیش‌بینی دوره های خشک و تر روزانه در کنیا استفاده‌کردند. آن‌هابا 
استفاده از الگوی زنجبره ما رکف مرتبه‌ی اوّل و توزیع احتمالاتی پ وآسون» طول دوره‌های خشکسالی و تر 
سالی بحرانی را در این منطقه تعیین کردند. 

عدل (۱۳۶۹) نیز از الگو زنجیره‌ی مار کف مرتبه‌ی اول» برای تعیین احتمالات تأمین آب از دریاچه‌ی 
سد امیر کبیر استفاده کرد. وی نتیجه گرفت که احتمال کمبود آب دریاچه در درازمدّت به کمک آمار 
۱ ساله» ۲./ است. 

در استان سیستان و بلوچستان به عنوان منطقه‌ی خشک و نیمه‌خشک ایران» بروز خشکسالی دور از 
انتظار نیست. این موضوع نشان‌دهنده‌ی آن است که برای انجام تحقیقات در زمینه‌ی خشکسالی باید تا کید 


1. Gabriel and Neumann 
2. Tel-Aviv 
3. Ochola and Kerkides 


سال هفتم پایش و پیش‌بینی خشکسالی در استان سیستان و بلوچستان ۳ 


بیشتری بشود تا نتایج به دست آمده آنها بتواند به عنوان ابزاری در برنامه ریزی‌های توسعه‌ی کشاورزی در 
این استان مورد استفاده قرا رگیرد. این که صبر شود تا خشکسالی به پایان برسد و سپس در مورد عواقب 
آن برنام‌ریزی کرد کاری غیر معقول و غیر علمی است و باید مدیریت ریسک را جایگزین مدیریت 
بحران نمود. در این مطالعه با استفاده از تحلیل الگویی که جهت پایش و پیش‌بینی خشکسالی تهیه شده 
است» وضعیّت خشکسالی در گذشته و آینده در استان سیستان و بلوچستان را مورد بررسی قرار می‌دهیم 


تا مدیریت‌ها و برنامه ریزی‌های لازم جهت کاهش ارات خشکسالی انجام شود. 


۲ مواد و روش‌ها 

استان سیستان و بلوچستان با وسعتی حدود ۱ کبلومتر مربع» ۱ درصد از مساحت کل کشور 
را تشکیل داده و از پهناورترین استان‌های کشور است و با قرار گرفتن در بین ۲۵ درجه و ۳دقيقه تا ۳۱ 
درجه و ۲۷ دقيقه عرض شمالی از خط استوا و ۵۸ درجه و ۵۰ دققه تا ۶۳ درجه و ۲۱ دقيقه طول شرقی از 
نصف‌النهار گرینویچ از نظر جمعیتی از کم ترا کم‌ترین استان‌های کشور است. این استان از شمال به استان 
خراسان جنوبی ‏ وکشور افغانستان» از شرق به کشورهای پا کستان و افغانستان از جنوب به دریای عمان و 
از مغرب به استان‌های کرمان و هرم زگان محدود می‌شود. بر اساس تقسیمات کشوری در سال ۸۱۳۸۱ 
استان سیستان و بلوچستان دارای ۳۶ بخش. ۳۱ شهر و ۹۸ دهستان می‌باشد(سالنامه آماری کشور ۱۳۸۱). 
استان سیستان‌وبلوچستان به لحاظ نزولات جوی اچیز و داشتن آب و هوای صحرایی از جمله استان های 
بسیار خشک کشور است. 


سیر 


شکل ۱. پرا کندگی مکانی ایستگاه‌های مورد مطالعه 


۴ مجّه‌ی جغرافیا و توسعه‌ی ناحه‌ای شماره‌ی یازدهم 

داده‌های اصلی این تحقیق را اطلاعات بارندگی ۸ایستگاه سینوپتیک استان سیستان و بلوچستان در 
طول دوره آماری ۳۰ ساله (۱۹۷۳-۲۰۰۳) تشکیل می‌دهد که از سازمان هواشناسی کشور تهیه شده‌است. 
مشخصات جغرافیایی این ایستگاه‌ها در جدول۱ ارائه شده است. نقشه‌ی منطقه‌ی مورد مطالعه و موقجّت 
ایستگاه‌ها در شکل ۱ آمده است. 

در نرم‌افزار تهیه o‏ با عنوان "DROUGHT MONITORING And FORECASTING"‏ که به 
اختصار ۴( نامیده شده‌است» برای درجه‌بندی خشکسالی از میان روش‌های موجود از دو شاخص !8۴ و 
درصد از نرمال که با توجه به داده‌های موجود در استان سیستان و بلوچستان کاربردی‌تر هستند» استفاده 
شده‌است وپیش‌بینی خشکسالی با استفاده از زنجیره‌ی ما رکف صورت گرفته‌است. در زیر به تشریح 


روش‌های مورد استفاده می‌پردازیم: 


جدول ۱. ایستگاه‌های سینوپتیکک مورد مطالعه و مشخصات جغرافیایی آنها 


ج | دت | در 
aT‏ 
a A KEN EA E‏ 


۱۳۳ r] س‎ 


مال‌هفتم...- پایش وپیش‌ینی خشکسالی دراستان سیستان وبلوچستان ها 
الف) شاخص ٩7‏ 

برای تعیین شاخحص استاندارد شده بارش از مقادیر بارند گی ماهانه‌ی هر ایستگاه در طول دوره آماری 
طولانی مدت (حداقل ۳۰ سال) استفاده شده‌است. ابتدا مقادیر بارند گی ماهانه‌ی هر ایستگاه برای هر یک از 
مقیاس‌های زمانی مورد نظر (...» ۵۲۴ ۱۲ ۶ ٩‏ ۳ ۱) محاسبه وسپس مقادیر بارندگی‌های تجمّعی در هر 
ماه به توزیع گاما برازش داده شده‌است. در نهایت این توزیع به یک توزیع نرمال تبدیل می‌شود. تابع 
توزیع گاما با استفاده ازتابع چگالی احتمالانی زیر محاسبه می‌شود (508062,2005): 


 -x/B‏ 0-1 ج 
BTC) e (۱)‏ ۳ = ))8 


در این رابطه: 
(20- تابع چگالی احتمالاتی 
6- مقادیر بارند گی تجمّعی در هر مقیاس زمانی و هر ماه از سال و برای هر ایستگاه هواشناسی مورد نظر (بر 
() 1 - تبع گاما 


0و 8 نیز پارامتر شکل و پارامتر مقیاس هستند که به صورت زیر محاسبه می‌شود: 


)۲( با )بت = 
/ اه 
a (۳‏ 
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که در آن: 
In‏ 
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در روابط بالا : 


× -میانگین بارند گی‌های تجمّی محاسبه شده در تمام ماههای مشابه در هر ایستگاه و در هر مقیاس زمانی 
-بارندگی تجمّعی در هر مقیاس زمانی و در هر ماه 


0- تعداد بارند گی‌های مخالف صفر ماههای مشابه در سری زمانی بارند گی‌های تجمّعی در هر مقیاس زمانی 


۶ مجل‌ی جغرافیا و توسعه‌ی ناحیه‌ای شماره‌ی بازدهم 
توزیع احتمال تجمّعی داده‌ها با انتگرالگیری روی توزیع گاما به دست می آید. 


۱ | ۵) 
01 )0( % 


G(x) = | ۵0۷ =‏ 
0 
6)0-احتمال تجمّعی توزیع گامای ناقص 
با توجه به این که توزیع احتمال تجمّعی گاما برای 0 = × تعریف نشده‌است و مسلماً در سری 
داده‌های بارندگی هر ایستگاه» مقادیر بارندگی صفر نیز وجود دارد لذا احتمال تجمّعی به صورت زیر 
تعریف می‌شود: 
H(x) = qt(1-qG() (۶)‏ 
٩-احتمال‏ وقوع بارند گی صفر در هر مقیاس زمانی» برای هر ماه و هر ایستگاه 
(116- توزیع احتمال تجمّعی گاما 
پس از محاسبه‌ی احتمال تجمّعی گاما در هر مقیاس زمانی و برای هر ماه از سال» این احتمال به بکك 
متغیر تصادفی نرمال استاندارد 2 با میانگین صفر و واریانس ۱ تبدیل می‌شود که این متفر تصادفی در 
حقیقت همان مقدار ٩۳1‏ مورد نظر است (1965, Abramowitz and Ste£Un‏ 


2 
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7,<5۳]- 4t ly, 0<H(K)<0.5 2 
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7 و6 


سال هفتم پایش و پیش‌بینی خشکسالی در استان سیستان و بلوچستان ۷ 


03<,ه 
c= 86‏ 
dı z.428‏ 
29 
58 -<:1 


مقدار 5۳1 مثبت نمایانگر بارش بیشتر از متوسط و مقدار 5۳1 منفی بیانگر بارش کمتر از مقدار متوسط 
است. زمانی که مقدار 5۴1 محاسبه‌شده منفی باشد. نشانه‌ی شروع خشکسالی است و هنگامی که این 
شاخص مثبت باشد» پايان خشکسالی را وید می‌دهد. شلات یک دوره‌ی خشکسالی؛ مجموع مقادیر ٩۳]‏ 
مثبت برای همه ماهها در همان دوره‌ی خشکسالی است. جهت تعیین درجه‌ی خشکسالی با استفاده از این 


شاعص از حدول۲ استفاده می شود: 


جدول ۲. مقادیر ٩۳1‏ و رطوبت و خشکی متناظر آن (2000 ,۲1۵765) 


مقاد بر 5۳1 


۸ مجلهی جغرافیاوتوسعه‌ی‌ناحیه‌ای ‏ شماره‌ی‌بازدهم ‏ 
ب) شاخص درصد از نرمال 

در تقسیم‌بندی شاخص درصد از نرمال ((۳) از توصیه‌های مهندسین مشاور کوانتا (۱۳۵۶) استفاده 
شده است. بر اساس این تقسیم‌بندی انواع شدّت‌های خشکسالی و ترسالی به صورت زیر ارائه گردیده 


است: 


PN = (P; / P)*100 ۱۱)‏ 
جدول۳. درجه بندی شلّت خشکسالی بر اساس روش درصد از نرمال( کوانتاه ۱۳۵۶) 
مقادبر درصد از نرمال(/) 
ترسالی شدید ۳۰< 
ترسالی خفیف ۱ ۳۰۱ 

نرمال 7-۰ ۱۰+ 

خشکسالی شدید ۱-۱ ۳۰- 
ی 5 ى ۰ فيو a‏ 

پ) زنجبره مار کف 


زنجیره‌های ما رکف در هیدرولوژی و هواشناسی جهت الگوسازی فرآیندهایی چون: بارندگی؛ 
جریان‌های رودخانه» رطوبت خاک و ذخیره‌ی آب در مخازن کاربرد دارد. اگر فرض شود که (26)0 
یک فرایند با مقدار گسسته باشد که در زمان صفر شروع شده و در طول زمان ادامه یابد» آن‌گاه مقادیر 


()×رابه صورت 2S‏ (...,0,1-) می‌توان نشان داد و داریم: 


۳12) = x [X(O) = xo, XO) = x... X(t-1) ی«‎ 9 
PXO =x [XCD =x] 


در رابطه‌ی(۱۲), ۲ احتمال شرطی است و بیان می کند که نتیجه‌ی هر فر آیند در زمان ا تنها به شرایط 
در زمان 1-]بستگی دارد. فرآیندهایی که چنین خاصیتی دارند یک زنجیره‌ی ما رکف مرتبه‌ی اول با 
زنجیره‌ی مار کف ساده نامیده می‌شوند (۲۳۵۳02500,1999), 

یک زنجیره‌ی مار کوف ساده با استفاده از ماتریس احتمال انتقالش ()۲ تعریف می‌شود. این ماتریس 


یک ماتریس مربع با عناصر ()/:9] است و به صورت زیر بیان می‌شود: 


۳] 2 -1( =i 2و‎ =j] 
p,)0 = ۳۵-۳ -(- 9 
P[X(t-1) =i] 


برای تمام جفت‌های آو ز, زنجیره ممکن است از حالت آدر زمان 1-1 به حالت ۲.. ,1,2 =[در 
زمان ] تغییر یابد, لذا با معلوم بودن وضعیّت در زمان 1-] احتمالات انتقال متناظر در زمان ]با 
2.۰۰ ورن نمایش داده می‌شود و لذا: 
i =1,2,..., (۴)‏ ۲,01 
j=1‏ 
چنانچه ماتریس احتمال انتقال (۳)0 وابسته به زمان نباشد» زنجیره‌ی مار کوف را ایستا یا همگن گویند 
و در این حالت می توان از علایم و P‏ به جای (۳)6 و () ص استفاده نمود )1993,ئMaidme1).‏ 


ماتریس احتمال انتقال به صورت زير نوشته می‌شود: 


2 | Po Po P22 Por (۱۵) 


1 Pro Pr مر‎ ۰ Pr ۱ 


۱۰ مجه‌ی جغرافیا و توسعه‌ی ناحیه‌ای شماره‌ی یازدهم 


عناصر قطر اصلی ماتریس بیانگر پایداری هر وضعیّت می‌باشد. 
چنانجه فرض شود زنجیره در حال حاضر در وضعیّت 1قرار دارد و پس از 10 مرحله‌ی زمانی در 
حالت زواقع شده است. احتمال گذر از حالت آبه حالت [در 1 مرحله را با 0 نشان داده و به صورت 


زیر بیان می‌شود: 


(n) : (n-1) ۳ نا‎ 
i Zk Fs ۶( 
Pij 2 Pk; ۳ 


چنانچه 0-1 آن گاه ر ,0 = )0 همان عناصر احتمال انتقال است. در رابطه‌ی (۱۶), 0/۳ عناصر 
ماتریس احتمال انتقال 10 مرحله‌ای ((۳۲۳) هستند که می‌تواند با ضرب کردن 1 بار ماتریس احتمال انتقال 
در خودش به دست آید. (2000 (Banik,‏ 


۳. روش تحقیق 

نرم افزار 01۴ برای محاسبات فوری پایش و پیش‌پینی خشکسالی تهیه شده‌است . ارائه‌ی این 
نرم‌افزار به دلیل گستردگی موضوع در مقاله ارائه نشده و علاقمندان برای دریافت آن می‌توانند با 
نگارند گان مکاتبه نمایند. با استفاده از نرم افزار 0۴ ۰ شاخص حساسیّت به خشکسالی ٩۳۱‏ 
درمقیاس‌های زمانی ۳ ۶ ۸ ۰۱۲ ۲۴ و ۴۸ ماهه برای تمام سال‌های آماری موجود محاسبه گردید و 
مقادیر سال ۲۰۰۳ آن برای مقباس‌های زمانی مختلف استخراج شد. با توجه به مقادیر ٩۳1‏ محاسبه شده 
برای سال‌های گذشته» درصد وقوع هریکک از انواع خشکسالی‌ها در مقیاس زمانی سه ماهه و برای هر 
یک از ایستگاه‌ها با توجه به فراوانی وقوع آنها در سی سال گذشته محاسبه شد (جدول۴. همچنین 
ماتریس‌های احتمال انتقال برای هشت ایستگاه استان سیستان و بلوچستان بر اساس شاخص 8۴1 (به 
صورت ماهانه) و شاخص درصد از ميانگین (به صورت سالانه) در یک دوره‌ی آماری سی ساله محاسبه 
شد. با مشخص بودن درجه‌ی خشکسالی یا ترسالی در هر ماه یا سال» می توان احتمال وقوع هر یک از 
شدّت‌های خشکسالی با ترسالی در ماه یا سال آتی را تعیین نمود (جداوّل۵الی ۸). این عمل با ضرب 
کردن مکرر ماتریس احتمال انتقال به دفعات مختلف در خودش به دست می آید که در نرم‌افزار 171۳ 


نیز وجود دارد. نقشه‌های پهنه‌بندی 5۲1 (برای ماه اول هر فصل.اشکال ۲الی ۴) و پیش‌بینی خشکسالی(به 


1. real time 


سال هفتم پایش و پیش‌بینی خشکسالی در استان سیستان و بلوچستان ۱۱ 
عنوان نمونه ماه ژانویه» اشکال ۵ الی ۸) برای منطقه‌ی مورد مطالعه با استفاده از نرم‌افزار 618و با استفاده از 


روش 1(۷ تهیه شده است. 


.٤‏ نتیجه گیری: 

باتوجه به نتایج به دست آمده از وضعیّت خشکسالی در سال‌های گذشته در منطقه‌ی مورد مطالعه 
چنین می توان نتیجه گیری نمود که» به طور کلی در یک مقیاس زمانی مشخص هر چه شدّت 
خشکسالی بیشتر می‌شود فراوانی حتمال وقوع آن نیز کاهش می‌یابد. در مقیاس‌های زمانی بلند مدّت 
(۱۲۰۹ و ...) به دلیل افزایش تعداد دفعات ضرب شدن ماتریس فراوانی درخودش درصد فراوانی انواع 
خشکسالی‌ها مقداری ثابت می‌شوده درواقع ماتریس تعادل منطقه به دست می‌آید. به طور کلی» منطقه‌ی 
مورد بررسی نسبت به خشکسالی‌های ملایم میان مدت و کوتاه مدت حساس‌تراست لذا می‌توان گفت 
که شاخص ٩۳‏ دارای یک روند منفی است و خشکسالی‌های اخیر باعث کاهش منابع آبی استان 
شده‌است. بنابراین شناخت همه جانبه‌ی حصوصیّات و ویژگی‌های خشکسالی ضرورت دارد تا بتوان به 
موفقیّت برنامه‌ریزی در مدیریت منابع آب کشاورزی امیدوار بود. با توجه به پیش‌بینی که با استفاده از 
زنجیره‌ی مار کف انجام شده است» مشخص شد که در ماه ژانویه سال ۲۰۰۴ احتمال وقوع خشکسالی 
کم در استان سیستان و بلوچستان به خصوص در مناطق م رکزی استان وجود دارد . احتمال وقوع 
خشکسالی در شهرستان‌های سراوان و خاش بیشتر است و مدیریت منابع آب باید با دقت بیشتری انجام 


گرفته و از توسعه‌ی منابع مصرف کننده در آن جل وگیری شود. 
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